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Abstract 

Recent results obtained by the KLOE experiment operating at DA$NE, the Frascati </>- 
factory, are presented. They mainly concern neutral kaon decays including the Kl dominant 
branching ratios, the Kl lifetime and the extraction of the CKM parameter V us from the 



Kl semileptonic decays and lifetime, 
is also presented. 



The best world upper limit on Ks — > ir n n channel 



1 Introduction 



The determination of \V US \ and \V u d\ provide the most precise test of CKM unitarity. In 
fact the first row unitarity requires \V u d\ 2 + \V US \ 2 + \V u b\ 2 = 1 which, since \V u b\ 2 ~ 10~ 5 , 
is equivalent to \V u d\ 2 + \ V US \ 2 = 1. The 2004 edition of Particle Data Group 1 gives \V u ,i\ = 
0.9738±0.0005 and \V US \ = 0.2200±0.0026 from which the sum of the squares gives 0.9966± 
0.0015 which deviates from unitarity by ~ 2a. Semileptonic kaon decays are the cleanest 
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way to obtain an accurate value of \V US \. In fact, since K — > 7r is a 0~ — > 0~ transition, only 
the vector part of the weak current has a non vanishing contribution and such processes are 
protected by the Ademollo-Gatto theorem against SU(3) breaking corrections to the lowest 
order in m s — ma- 

In order to extract \ V US \ from semileptonic kaon decays we need to measure the branching 
fraction and the Kl lifetime. The most accurate determination of \ V US | from Kl semileptonic 
decays comes from KTeV collaboration 2 . The K e z branching fraction is found to be 0.4067± 
0.0011, strikingly different from PDG fit value, 0.3881 ± 0.0027. The KTeV result has been 
partially confirmed by the NA48 Collaboration 3 , 0.4010 ± 0.0045, even if with an error <~ 4 
times bigger. 

The Kl lifetime value in the PDG relies on a single direct measurement tl = (51.54±0.44) 
ns that was performed more than 30 years ago 4 . Kl lifetime is, at present, the major 
experimental source of error in the determination of \V US \. 

KLOE has the unique possibility to measure simultaneously Kl absolute branching frac- 
tions and lifetime. The measurement of the Kl absolute branching is possible at the ^-factory 
since the production of KsKl pairs in <j> decays provides a tagged, monochromatic Kl beam 
of known flux. Moreover Kl have low momentum ~ 110 MeV/c) and, therefore, a 

big fraction (<~ 50%) of them decays inside the detector. 

This is important since the statistical error on the lifetime depends on the number of 
events, cx l/y/N and, very strongly, on the time interval covered: 



V (stat) = T = 

where T = St/r is the interval covered by the fit in lifetime units and N is the number of 
events in that interval. With T <~ 0.4 and N ~ 9 • 10 6 (which is the KLOE case) a statistical 
error of ~ 0.3% can, in principle, be reached. 
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2 Experimental Setup 

DA$NE, the Frascati <j> factory, is an e + e~ collider working at W ~ <~ 1.02 GeV with 
a design luminosity of 5 x 10 32 cm~ 2 s _1 . The <j) mesons are produced, almost at rest, with 
a visible cross section of ~ 3.2 ^b and decay into K + K~ (KsKl) pairs with BR of ~ 49% 
(<~ 34%). These pairs are produced in a pure J PC = 1~~ quantum state, so that observation 
of a Ks (K + ) in an event signals (tags) the presence of a Kl (K~) and vice- versa; highly 
pure and nearly monochromatic Ks, Kl,K + and K ~ beams can thus be obtained. Neutral 
kaons get a momentum of <~ 110 MeV/c which translates in a slow speed, (3k ~ 0.22. Ks 
and Kl can therefore be distinguished by their mean decay lengths: As ~ 0.6 cm and A^ <~ 
340 cm. 

The KLOE detector consists of a drift chamber, DCH, surrounded by an electromagnetic 
calorimeter, EMC. The DCH 5 is a cylinder of 4 m diameter and 3.3 m in length which 
constitutes a large fiducial volume for Kl decays (1/2 A^). The momentum resolution for 
tracks at large polar angle is a p /p < 0.4%. The EMC is a lead-scintillating fiber calorimeter 6 
consisting of a barrel and two endcaps which cover 98% of the solid angle. The energy resolu- 
tion is oe/E ~ 5.7%/ v /E(GcV). The intrinsic time resolution is a T = 54 ps/^/E(GeV)® 50 
ps. A super-conducting coil surrounding the barrel provides a 0.52 T magnetic field. 

During 2002 data taking, the maximum luminosity reached by DA$NE was 7.5xl0 31 
cm~ 2 s _1 . Although this is lower than the design value, the performance of the machine was 
improving during the years and, at the end of 2002, we collected ~ 4.5 pb _1 /day. The whole 
data sample in the years 2001-2002 amounts to 450 pb _1 , equivalent to 1.4 billion <j> decays. 
The analyses presented here are based on ~ 400 pb _1 of integrated luminosity of the 2001 
and 2002 runs. 

Recently, the machine has been upgraded and KLOE has resumed its data taking in 
April 2004. Up to know (15 th May, 2005) <~ 1 fb _1 have already been collected with a peak 
luminosity of 1.3xl0 32 cm -2 s _1 . We foresee to reach ~ 2 fb _1 by the end of the year. 



3 Measurement of the dominant Kl branching ratios 



The Kl absolute branching fractions are determined on a tagged Kl events sample by 
counting the number of Kl decays in each channel, A^, in the used fiducial volume and 
correcting for acceptance, reconstruction efficiency and tagging efficiency: 

N- 1 
BR(K L ^i) = —^x — — — (2) 

Ntag frecO) X e F V\T L ) X 6 tag (l) / C tag (all ) 

where i = 7r°7r 7r°, ir e^v, ir [i^v, tt + ir~ 7r° , N tag is the number of tagging events, e rec is the 
reconstruction efficiency (~ 45% for ir + ir~ir° events, ~ 60% for n e^v, n fi^u events and 
~ 100% for 7r 7r°7r° events), €fv(tl) is the fiducial volume geometrical acceptance which 
depends on the Kl lifetime and e ta g(i)/etag(all) is the fractional variation of the tagging 
efficiency when the Kl decays in a given channel with respect to the average. We define this 
ratio tag bias. 



3.1 The tag 

The tag is provided by K$ — * ir + n~ events selected requiring the presence of a vertex with 
two opposite curvature tracks within a cylinder of radius r < 10 cm and height h < 20 cm 
around the interaction point (IP). The two-tracks invariant mass, in the pions hypothesis, 
must be within 5 MeV around mx s ■ The magnitude of the total momentum of the two tracks 
must be within 10 MeV of the value expected from the value of p^. The Kl momentum 
is obtained from the decay kinematics of (f> — > KsKl using the Kg direction reconstructed 
from the measured momenta of ir + ir~ tracks and the known value of p^. 

The main source of tag bias is due to the dependence of the trigger efficiency on the Kl 
behaviour. The hardware calorimeter trigger which requires two local energy deposits above 
some thresholds (50 MeV on the barrel and 150 MeV on the end caps) is used for the present 
analysis. The trigger efficiency is essentially 100% for 7r 7r°7r°, between 95 — 85% for charged 
decays. To reduce the tag bias due to trigger efficiency we require that the trigger conditions 
are satisfied by the pions from K$- We further reinforce these conditions by requiring that 
the two pions impinge on the calorimeter barrel and produce two clusters with an energy 
E > 80 MeV each. 

The FV used for the analysis is defined inside the drift chamber by 35 cm < \J x 2 + y 2 < 
150 cm and \z\ < 120 cm, where (x, y, z) are the Kl decay vertex position coordinates. Since 
the Kl mean decay path length in KLOE is <~ 340 cm, the FV contains <~ 26.1% of the 
Kl decays. This choice minimize the difference in tag bias among the decay modes. The 
average tag bias is 0.985, 0.99 and 1.02 for 7r ± e =F ^ or 7r ± ^ =F ^, 7r + 7r _ 7r° and 7r 7r°7r° decays 
respectively. 



3.2 Charged Kl decays 

The Kl in charged decay modes are selected by requiring the presence of two good tracks 
forming a vertex in the FV and not belonging to the Kg decay tree. A track is associated to 
the Kl if the point of closest approach to the Kl line of flight has a distance with respect 
to the Kl line of flight d c < a^J x 2 + y 2 + b, with a — 0.03 and b— 5 cm. 

The tracking efficiency has been determined by counting the number of events with at 
least one found Kl track and the number of events in which there are two opposite sign 
decay tracks. The tracking efficiency is also evaluated from Monte Carlo simulation and 
it is 60.5% for K e3 , 58.5% for K^ and 43.0% for 7r + 7r~7r°. A correction to the tracking 
efficiency has been applied by comparing the results for data and Monte Carlo simulation. 
The correction is evaluated as a function of the track momentum using Kl — * 7r + 7r~7r and 
Kl — > 7r ± e T ^ events and ranges between 1.03 and 0.99 depending on the channel. 

The variable to discriminate among the different charged decay modes is the lesser of the 
two values of = |p mis | - E mlss , where 

Pmis is the missing momentum and E m i SS is the 
missing energy evaluated in the two mass assignments 7r + , /i - or 7r~/i + . An example of this 
distribution is shown in Fig. 1 where the different components are shown. 




Figure 1: A TM distribution for an event subsample. Dots are data, solid histograms are Monte Carlo. 




The A,,-^ distribution obtained with data is fitted with a linear combination of three 
Monte Carlo distributions (for it e^v, n^/x^v and 7r + 7r~7r° events) by leaving free the relative 
weights. The contribution from the CP violating decay Kl — > 7r + 7r~ and the Kl ~ > 7r 7r 7r° 
with Dalitz conversion is kept fixed in the fit and amounts to 0.3 %. In Fig. 2 we show 
enriched samples of K e % (left) and K^ (right) events obtained by selecting the e or /i by 
time of arrival and energy deposition in the calorimeter. The radiative corrections which 
affect mainly the K e % decays have been properly taken into account in the Monte Carlo 
simulation 7 . 

3.3 Neutral Kl decays 

Kl — > 7r 7r°7r° decays are selected using time of flight techniques. In fact, the position of 
the Kl vertex for Kl — > 7r 7r°7r° decays is measured using the photon arrival times on the 
EMC. Each photon defines a time of flight triangle shown in Fig. 3. The three sides are the 
Kl decay length, Lk; the distance from the decay vertex to the calorimeter cluster centroid, 
L 7 ; and the distance from the cluster to the <f> vertex, L. The equations to determine the 
unknowns Lk and L 1 are: 
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Figure 3: Left: the time of flight triangle. 



L 2 + L K -2LL K cos8 = L 2 (3) 
Lk 1 0K + Ly = ci 7 

where i 7 is the photon arrival time on the EMC, @k c is the Kl velocity and 9 is the angle 
between L and Lk- Only one of the two solutions is kinematically correct. The Kl vertex 
position is obtained by the energy weighted average of each Lk measurement. 

The accuracy of this method is checked with Kl — > 7r + 7r~7r° decays, comparing the 
position of the Kl decay vertex from tracking using the 7r+7r~ pair with the one from timing 
with the two photons from tt°. The vertex reconstructed by the calorimeter has on average 
an offset of 2 mm almost uniform in the fiducial volume. 

To select the Kl — * 7r 7r°7r° events we require at least three photons with energy greater 
than 20 MeV originating from the same vertex. 

The main sources of inefficiencies are: 1) geometrical acceptance; 2) cluster energy thresh- 
old; 3) merging of clusters; 4) accidental association to a charged track; 5) Dalitz decay of 
one or more 7r°'s. The effect of these inefficiencies is to modify the relative population for 
events with 3, 4, 5, 6, 7 and > 8, clusters with a loss of efficiency of ~ 0.8%. 

Background contamination affects only events with three and four clusters. The main 
source of background comes from Kl — ► 7r+7r~7r° decays where one or two charged pions 
produce a cluster not associated to a track and neither track is associated to the Kl vertex. 
This background is rejected by requiring at least one cluster with at least 50 MeV energy 
and a polar angle satisfying |cos#| < 0.88. Other sources of background are Kl — > ir°n 
decays, possibly in coincidence with machine background particles (e ± or 7) that shower 
in the QCAL and generate soft neutral particles and K$ — > 7r°7r° following the Kl — > Ks 
regeneration in the drift chamber material. 

34 Results 

A total of ~ 13 milions of tagged Kl are used to compute the branching fractions, almost ~ 
40 milions to evaluate systematic uncertainties. Using the published result of the Kl lifetime 
(jl = 51.54 ± 0.44 nsf we obtain the following results: 
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Table 1: Summary of systematic uncertainties on the absolute branching fractions measurements. 



0.4049 ± 0.0010 stat ± 0.0031 syst (4) 
0.2726 ± 0.0008 stat ± 0.0022 syst 
0.2018 ± 0.0004 stat ± 0.0026 syst 
0.1276 ±0.0006 stat ±0.0016 syst 

(5) 

where the sources of systematic uncertainties are shown in Table 1. The sum of all measured 
branching fraction above, plus the PDG value for rare decays, 0.0036, is ^BRi = 1.0104 ± 
0.0018 corre i a |- ec | ± 0.0074 uncorre i a -j- e( j where the correlated error includes all contributions 
to the uncertainties on the branching ratios that are 100% correlated between channels, such 
as the uncertainty in the value of the Kl lifetime. This result depends on the value of the 
Kl lifetime through the acceptance (Eq. 2). Turning the argument around, by normalizing 
the sum to 1 we obtain un indirect estimate of the Kl lifetime: 

t(K l ) = (50.72 ± 0.14 stat ± 0.36 syst )ns 
and a new set of values for the branching fractions: 



BR(K L -» n ± e Zf: i>) = 

BR(K L -> tt^i/) = 

BR(K L -» ttW) = 

BR(K L -» TT+TT-TT ) = 



BR(K L -> n ± e T v) = 0.4007 ± 0.0006 stat ± 0.0014 syst (6) 
BR(K L -> tt ± ij^i/) = 0.2698 ± 0.0006 stat ± 0.0014 syst 
BR(K L -» 7r 7r°7r°) = 0.1997 ± 0.0005 stat ± 0.0019 syst 
BR(K L -» tt+tt^tt ) = 0.1263 ± 0.0005 stat ± 0.001 l sy;st 

(7) 



4 Direct measurement of the Kl lifetime 

We have measured the Kl lifetime using <~ 15 x 10 6 events of the fully neutral decay Kl — » 
7r 7r°7r° tagged by Kg — » 7r + 7r~ events. This choice is motivated by the fact that we want 
to maximize the number of tagged events to reduce the statistical error and, simultaneously, 
we want to avoid any coupling among tagging and tagged events in order to minimize the 
systematic uncertainty. 

The Kl — > tt°tt tt° sample selected for the BR measurement has been used also for the 
Kl lifetime measurement plus some additional cuts. 

For lifetime measurement we must keep under control the variation of the selection ef- 
ficiency with Lk- Monte Carlo simulation shows that the selection efficiency has a linear 
dependence with Lk, e{L K ) = (0.9921 ±0.002) - (1.9±0.2) • 10" 5 -^(cm) mainly due to the 
vertex reconstruction efficiency. The vertex reconstruction efficiency as a function of Lk has 
been checked also using Kl — * tt + it~it events both in data and in Monte Carlo simulation. 
We find the same linear dependence as in the Kl — * 7r 7r°7r° case with slopes compatibles 
within their statistical uncertainties. 
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Table 2: Left: fit of the proper time distribution. Right: residuals of the fit. 



The Kl lifetime is measured using events with a vertex reconstructed in the region 40 
cm < Lk < 165 cm and with a flight direction defined by a polar angle 8 with respect to 
the beam axis between 40° < 9 < 140°. These two conditions define the fiducial volume. 

The Kl proper time, t* , is obtained event by event dividing the decay length Lk by /?7 
of the Kl in the laboratory, t* = Lfc/(/3'yc). The residual background is subtracted bin by 
bin using Monte Carlo predictions. 

The variations of the vertex reconstruction efficiency as a function of the decay length are 
taken into account by correcting bin by bin the decay vertex distribution with the efficiency 
values obtained with the Kl — > 7r 7r 7r° Monte Carlo sample multiplied by the data - Monte 
Carlo ratio of the efficiencies evaluated with 7r + 7r _ 7r° sample. 

The statistical uncertainty of the efficiency values (~ 0.1%) has been taken into account 
by adding it in quadrature to the statistical fluctuation of the entries in each bin of the t* 
distribution after the background subtraction. 

The distribution is fitted with an exponential function inside the fiducial volume, which, 
in terms of proper time ranges from 6 ns to 24.8 ns. This corresponds to a time interval 
T ~ 0.37 expressed in lifetime units. With ~ 8.5 • 10 6 events inside the fit region we obtain 
t = (50.87 ± 0.17 s t a t) ns with a % 2 = 58 for 62 degrees of freedom. The major sources 
of systematic uncertainties come from the background evaluation, tagging and selection 
efficiency and from the estimate of the Kl nuclear interactions in the drift chamber material. 
The total systematic error is ~ 0.5%. Our result is: 

t(K l ) = (50.87 ± 0.17 stat ± 0.25 syst )ns 

ft is compatible at 1.3 a level with the other measurement 4 and only at 1.7cr level with the 
PDG 2004 fit 1 . 



5 Determination of \V U 



\V US \ is proportional to the square root of the semileptonic BR of K mesons. For Kf^ decays 



we can write 



8. 
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where f+°(Q) is the vector form factor at zero momentum transfer and Ii(X + ,X' + ) is the 

integral of the phase space density, factoring out f^"'" and without radiative corrections. 
Radiative corrections at large scale of form factor and phase space density are contained in 



the term = (0.55 ±0.10)% 9 ' 10 . 



The short-distance electroweak corrections are included 



A + and X', are the slope and curvature of the vector form factor 



in the parameter S e 

/£ 71 (0). From eq. 8 we see that the lifetime value enters directly in the determination of 
the product \V US \ x f¥ O7r (0). The three most recent results are: 



BR{K L 
BR(K L 
BR(K L 



r±P=F, 



= (40.67 ±0.11)% KTeV 
= (40.10 ±0.45)% NA48 
= (40.07 ±0.15)% KLOE 



where the error is the sum in quadrature of the statistical and systematic uncertainties. Fig. 
4 (left) shows the product \V US \ x w (0) extracted from the three recent measurements of 
BR(K L -> i^e^v) using the PDG 2004 average r = (51.50 ± 0.40) ns. 

In the extraction of V us we use the values of A + and X' + obtained by KTeV experiment 
from a quadratic fit 11 . In the same plot we show also the value of \V US \ x f+ OjT (0) obtained 
from the KLOE preliminary measurement 12 of BR(K$ — > if^e^v) = (7.09 ± 0.07 s t a t ± 
0.08 syst ) ■ 10~ 4 and the ±ltr band from the \V u d\ and unitarity where we use f? "(O) = 
0.961 ± 0.008 following Ref. 13 . Fig. 4 (right) shows the same product extracted using the 
KLOE K L lifetime value, r L = (50.87 ± 0.17 ± 0.25) ns. The K L data are now in better 
agreement with K$ ones and unitarity. 
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Figure 4: \V US \ x fF "'(O) for various K semileptonic decays. 



6 Direct search of Ks — > 7r 7r°7r° decay 

The decay K$ — > 37r° is a pure CP violating process. The related CP violation parameter 
77000 is defined as the ratio of K$ to Kl decay amplitudes: 77000 = A{Ks — + 3-k°)/A{Kl — > 
£ 000 where e describes the CP violation in the mixing matrix and e 000 is a direct 
CP violating term. In the standard model we expect 77000 to be similar to 7700 (|??oo| ~ 
2 x 10~ 3 ). The expected branching ratio of this decay is therefore ~ 2x 10~ 9 , making its 
direct observation really challenging. The best upper limit on the BR (i.e. on | ?7ooo 1 2 ) has 
been set to 1.4xl0 -5 by SND 14 where, similarly to KLOE, it is possible to tag a K$ beam. 

The other existing technique is to detect the interference term between KgK^ in the 
same final state which is proportional to 77000 , Sfyooo cos(Amt) — ^77000 sin(Amt). The 
best published result using this method comes from the NA48 Collaboration 15 . Fitting 
the Ks — Kl — * 7r°7r°7r interference pattern at small decay times, they find 3?77ooo = 
-0.002 ± O.Ollstot ± 0.015 syst and ^77000 = -0.003 ± 0.013 stQ t ± 0.017 syst corresponding 
to BR(K S -> 37r°7r 7r°) < 7.4 x 10~ 7 at 90% C.L. 

The signal selection requires a Kl interacting with the calorimeter and six neutral clusters 
coming from the interaction point (IP). A first rejection of the huge background coming 
from the decay Ks — > ir°ir + 2 fake 7 is obtained by applying a kinematic fit imposing as 



constraints the K$ mass, the Kl 4-momcntum and (3 = 1 for each photon. Two pseudo-^ 2 
variables are then built: x^o which is based on the best 67 combination into 7r 7r°7r° and 
X% n o which selects four out of six photons providing the best kinematic agreement with the 
Ks — > 7r°7r° decay. A signal box is defined in the x\^o vs x^o plane by optimising the upper 
limit in the Monte Carlo sample. Residual background comes from Ks — > tt + it~ , Kl — > 37r° 
events in which one of the two pions interacting with the quadrupoles produces a late cluster 
simulating a AT^-crash. This background is rejected by vetoing events with two charged 
tracks from IP. 

At the analysis end we find 2 events in the signal box with an estimated background of 
B = 3.13 ± 0.82 sta t ± 0.37 sy st . To derive the upper limit on the number of signal counts, 
we build the background probability distribution function, taking into account our finite 
MC statistics and the uncertainties on the MC calibration factors. This function is folded 
with a Gaussian of width equivalent to the entire systematic uncertainty on the background. 
Using the Neyman construction 16 we limit the number of Ks — > 7r°7r 7r° decays observed 
to 3.45 at 90% C. L. with a total reconstruction efficiency of (24.36 ± 0.11 sta t ± 0.57 syst )%. 
In the same tagged sample, we count 3.78 x 10 7 Ks — > tt°tt° events. We use them as 
normalization. Finally, using the value BR(Ks — > 7r°7r°) = 0.3105 ± 0.0014 1 we obtain 
BR(Ks — > 7r 7r°7r°) < 1.2 x 10~ 7 at 90 % C.L. which represents an improvement by a factor 
~ 6 with respect to the best previous limit'. 
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